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ABSTRACT 

Recently the Milagro experiment observed diffuse multi-TeV gamma ray emis- 
sion in the Cygnus region, which is significantly stronger than what predicted by 
the Galactic cosmic ray model. However, the sub-GeV observation by EGRET 
shows no excess to the prediction based on the same model. This TeV excess 
implies possible high energy cosmic rays populated in the region with harder 
spectrum than that observed on the Earth. In the work we studied this theo- 
retical speculation in detail. We find that, a diffuse proton source with power 
index ctp < 2.3, or a diffuse electron source with power index Ue ^ 2.6 can re- 
produce the Milagro's observation without conflicting with the EGRET data. 
Further detections on neutrinos, a diagnostic of the hadronic model, and hard 
X-ray synchrontron radiation, a diagnostic of the lepton model, help to break this 
degeneracy. In combination with the gamma ray observations to several hundred 
GeV by Fermi, we will be able to understand the diffuse emission mechanisms in 
the Cygnus region better. 

Subject headings: ISM: general — cosmic rays — acceleration of particles — 
gamma rays: theory 



Introduction 



The Galactic diffuse gamma-ray emission provid es important information on the origin 
and propagation of the Galactic cosmic rays (GCRs) (IStrong et al.ll2000l . l2004bl ). In general, 
there exist three possible mechanisms for the diffuse gamma-ray emission: decay of neutral 
pion produced by hadronic interaction between cosmic ray (CR) nuclei and interstellar gas; 
inverse Compton (IC) scattering between CR electron and interstellar radiation field (ISRF); 
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and the bremsstrahlung radiation by interaction between CR electron and interstellar gas. 
In addi tion, dark matter annihilation is eraerging as an alternative possi bility to the diffuse 
7-rays Jde Boeilboosi : Ide Boer et aUboOsI : Ide Boei][2007l : IbI et aljl2008al lb[). 



The Cygnus region, defined as 65° < / < 85° and —3° < 6 < 3° following the Mi- 
lagro measurements, is rich in molecular clouds and is one of the richest star formation 
regi ons in the Galaxy (|Dobashi et al.l Il996l: ISchneider et al.M2006r). Observations in ra- 
dio (I 



Piddington &: Minnett 



20031 ). optical (IDickel et al. 



1952), infrare d (IKnodlsedei 



19691 ). X-ray (iGiacconi et al. 



20o3 iComeron et al 



19671 ) and 7-ray (IChen et al. 



2002; 



Hanson 



19961 ) 



bands found many interesting sources in this region. These sources have provided us valu- 
able information about the astrophysical processes in this region. Ground-based very high 
energy (VHE) 7-ray observatories also det ected TeV 7-ray emissions in the Cygnus region 
( lAharonian et al.l l2002l : lAbdo et al.ll2007a] ). Such detection is of great importance for CR 
physics, since it indicates the existence of high energy CR accelerators. Especially the Mila- 
gro experiment^ has made remarkable progress recently. 



(Abdo et al. 


2007a. 


2008. 


2007b) 



J2031+41 is possibly th e counterpart of the un identified TeV source TeV J2032+4130, first 



discovered by HEGRA flAharonian et al.l 120021 ) . and MGRO J2019-F37 is observed for the 



first time without an obvio us counterpart. There are extensive discussions about the na- 



ture of the two TeV sou rces (lAharonian et al.ll2002l : iButt et al.ll2003l : iMukherjee et al.ll2003 



Beacom fc Kistlerll2007l ). However, no consensus has been reached yet. As for the diffuse 
7-ray emission, the Milagro measurement shows an evident excess (hereafter "TeV excess") 
compared with the diffuse emission predicted by the Galactic cosmic ray model GALPROP 
(~ 7 tir nes of the convent ional model prediction and ~ 2 times of the optimized model pre- 
diction, |Abdo_et_alJ|2008|). On the other hand, the s ub-GeV band measurements by EGRET 
are well consistent with the GALPROP prediction (IHunter et al.l Il997l : lAbdo et al.l l2007al ) 

B 

The "TeV excess" may be due to some unident ified TeV sources o r diffuse high energy 
proton or electron population in the Cygnus region (lAbdo et al.ll2007al ). If the latter is true 
we can draw two general conclusions according to the Milagro and EGRET data: 1), the 



^Milagro homepage, |http://umdgrb.uind.edu/cosmic/milagro.html 



■^It should be pointed out that, there are "GeV excesses" of E GRET observations compared with the 
conventional model of diffuse 7-ray emission in >GeV energy band (jHunter et al.lll997r ). However, since it 
is not the main purpose of the present work, we just briefly discuss the "GeV excess" problem for the sake 
of completeness and do not intend to go deep in this topic. 
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spectrum of the proton or electron should be harder than the one observed at the Earth; 
2), the CRs of the source population need to reside in the Cygnus region in order not to 
exceed the locally measured CR fluxes. In this work we investigate this possible explanation 
in detail by building realistic models and deriving their implications for future experiments. 

This paper is organized as follows: in Sec. 2 we give an introduction to the diffuse 7-ray 
emission predicted by conventional GCR model. In Sec. 3 we present a brief introduction 
to the "GeV excess" problem. Our models to solve the "TeV excess" problem and possible 
observational effects for future experiments are discussed in Sec. 4. Finally the conclusion is 
drawn in Sec. 5. 



2. Diffuse 7-rays from the Conventional GCR model 



The GCRs propagate diffusively in the Galactic magnetic field (see e.g.. lGaisserlll990l ). 
The interactions between CRs and the interstellar medium (ISM) can generate diffuse 
rays during propagat i on. A numerical solution of the diffusion equation is developed by 
Strong fc Moskalenkol JlQQsh . which is known as the GALPROP model. In GALPROP the 
model parameters are adjusted to reproduce the locally measured CR spectra. The realistic 
distributions of ISM and ISRF are incorporated in GALPROP to calculate the fragmenta- 
tions and energy losses. Since the realistic astrophysical inputs and the general treatments 
of relevant physical processes, such as convection and reacceleration, the GALPROP is rec- 
ognized as the best CR propagation model at present. 

We adopt the conventional GALPROP model to calculate the GCRs propagation and 
the 7-ray emission (hereafter we denoted this component as "GC R background"). The 



isotropic extra Gala ctic (EG) component of diffuse 7-rays is included (ISreekumar et al.lll998 



Strong et al.ll2004al ). The results are shown in Fig HI The EGRET data are adopted from 
the d iffuse sky maps with the point sources from 3EG catalog subtracted (jCillis fc Hartman 
2OO5I ). The figure shows that the GALPROP prediction matches well with the low energy 
measurements of EGRET data (except the "GeV excess" ), while for the TeV data of Milagro, 
it shows a remarkable under-estimation. The model prediction is about 6 times lower than 
observation. This is consistent with the results given in lAbdo et al.l (120081 ). It has been shown 
that even the optimized GALPROP model, which is developed to solve the E GRET "GeV 



excess" prob lem (see the next section), predicts only a half of the Milagro data (lAbdo et al. 



2007al . l2008l ). Therefore there must be an additional photon component with harder spectrum 
in the Cygnus region. We also note that the local high density of ISM or ISRF in the Cygnus 
region can not account for this "TeV excess" without violating the EGRET observations, 
since it just makes the curves shift upwards. 
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Fig. 1. — The diffuse 7-ray spectrum in the Cygnus region predicted b y the conventional 
GAL PROP model compared with the observational data from EGRET ( ICillis fc Hartman 



20051) and Milagro (lAbdo et al.ll2007ah . 



3. The GeV excess 



The Galactic diffuse 7-rays measured by EGRET show an excess for energies ^GeV 
compared with p redictions by the c onventional GCR model, which is referred as "GeV 
excess" problem (IHunter et al.l 119971 ) . The "GeV excess" appears in all directions of the 
sky, and there is no special structure in the Cygnus region. Many models are proposed 
to solve the "Ge V excess" problem by t uning the C R spectra, such as the harder nucleon 



spectrum model f|Gralewicz et al. 



( Strong et al 



( IStrong et al 



2000[). However, these models usually cannot reproduce the all-sky 7-ray data 



1997 



Mori! 119971 ) ■ the harder electron spectrum model 



( Porter fc Protheroel 1997 : Pohl fc Espositol " 19981 ) or the combination of the two models 



2004a). 



Strong et al.l (j2004bl ) have developed an "optimized" model by adjusting the interstellar 



proton and electron intensities to reproduce the 7-ray data and CR data, such as B/C, 
simultaneously. The optimized model gives a good fitting to the EGRET data of all sky 
directions. However, large spatial fluctuations of proton and electron fluxes have to be 
incorporated in their model (jStrong et al.ll2004bl ). 
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Another approach to the "G eV excess" problem is by DM annihilation (jde Boerl 12005 



de Boer et al. 


2005; 


de Boer 


2007) 



with predictions after taking the DM annihilation into acc o unt if assuming; a supersymmetric 
DM with mass m^^ - 50 - 70 GeV Jde Boer et aljl2005h . IbI et all J2008al lb[) extended the 
DM scenario using DM subhaloes to account for the "boost factor" and calculate both the 
CR background and D M signals in mor e realistic propagation models. In the present work 
we adopt the model of iBi et al.l (l2008al ) to fit the GeV data. However, it should be noted 
that it is far away from the final answer to the "GeV excess" problem. To go details of this 
issue is beyond the scope of the present study. 



TeV excess 



In this section we turn to the "TeV excess" problem. We introduce a high energy 
population of CRs to explain the Milagro data. It is known t hat the Cygnus reg ion is 

booik OB 



rich in poten tial CR accelerators, such as the Wolf-Rayet stars (Ivan der Huch 
associations ( Bochkarev &: SitniklllQSsI ) and supernova remnants (SNRs, 



Green 



20041 ). Here 



we assume there exist a (or several) CR accelerating source (s) in this region. The ages of the 
sources should not be too old (e.g., 10^ ~ 10^ yr, for typical diffuse coefficient D ~ 10^^ cm^ 
s~^, the propagation length after this time is about several hundred pc) so that the CRs can 
distribute diffusively in the Cygnus region and will not affect the measurements of CR fluxes 
at the Earth. We discuss in detail the hadronic and leptonic scenarios in the following. 



4.1. The hadronic model 



We flrst consider a proton source uniformly distributed in the Cygnus region. We assume 
it has the energy spectrum as d0p/di?p oc Ep"^ exp{—Ep/Ep). The spectrum index Up is 
adopted as 2.3, which can well reproduce the Milagro data and is still consistent with the 
low energy data by EGRET. It should be noted that the adoption of Up is not arbitrary. On 
one hand, it cannot be too soft in order not to exceed the sub-GeV observations by EGRET. 
On t he other hand, the Fermi acceleration of CRs from shock waves predicts spectrum index 
~ 2 (jGaisseiill990l ): and the diffusion in the interstellar medium can soften the spectrum a 
bit (e.g., 1/3 for a Kolmogorov diffusion). Therefore the Milagro observations can indeed 
limit the source spectrum in a relative narrow range. The cut-off energy Ep is due to the 
acceleration limit of the source. For supernova blast shock it is known to be 100 ~ 1000 
TeV for protons (IBerezhkol Il996l ). The shock acceleration by stellar winds of red giants 
and Wolf-Rayet stars gives 10 PeV (IVoelk &: Biermannl 119881 ). The ensemble shocks of OB 
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association can even give the maximum energy about 100 PeV (IBykov fc ToptyginI Il990l ). 
Here we adopt the cut-off energy to be 1 ~ 10 PeV. Such a cut-off seems not to affect the 
results of 7-rays at ~ 10 TeV. However, it will affect the high energy neutrino ffux, as shown 
below. The normalization of the CR proton density is tuned to be consistent with the high 
energy observations by Milagro. 

The 7-ray spectrum after taking such new sources into account is shown in Figj2l It 
shows a good agreement between the model prediction and observations. We notice that 
the high energy behavior of 7-rays is dominated by the new proton source contribution, so 
precise measurements of the high energy spectrum (from several GeV to hundreds of TeV) 
will be helpful to determine the spectrum of the source term. The total energy of the proton 

)( i kpc )^ ^^S- For typical ISM density 
1 kpc (which is corn parable with the 



51 / 1 cm 



source is estimated as W{> 1 GeV) ~ 10 ( 
uh ~ 1 cm^"^ and the distance of the source d 
molecular clouds and OB associations in this region, [Schneider et al.ll2006r). the tota l energy 
is similar with the energy release of a typical core-collapse supernova (jZwickylll962l ). 




Fig. 2. — The energy spectrum of diffuse 7-rays in the Cygnus region as measured by EGRET 
and Milagro and predicted by the GALPROP model combined with additional proton source 
and DM annihilation components. See text for details. 



For inelast ic p — p collis ion, the pion mesons 7r+, tt and 7r° are generated with almost 
equal amount (IGaisserl Il990h . 7-rays are produced through the decay of 7r° mesons. The 



- 7- 



decay of tt^ will lead to neutrino emission. In general the initial neutrino flavor ratio is 



becomes : v., 



2:0. Because of neutrino oscillation the flavor ratio on the Earth 
= 1:1:1. The ratio between the total number of neutrinos and 7-ray 



photons is 3 : 1. The t ypical energy of neutrinos is 



flBeacom fc Kistleil 120071 ) 



1/2 of the 7 photons, so we have 



dE, 



, #7 



for each kind of neutrinos. 



In Fig. [3] we show the cumulative muon flux induced by + on a km^ detector, 
like IceCube, for one-year observation. The muon flux including both the containing events 
(that muons generated in the detector) and the through-going events (that muons generated 
outside th e detector and pro p agate s into the detector volume) is calculated following the 
method of iKistler fc BeacomI (120061 ). The absorption of neutrinos by the Earth along the 
Cygnus direction (R.A. ~ 20^*, Dec. ~ 40°) is taken into account. Thi s absorption effect 
is known to be significant for neutrino energies greater than ~ lOTeV (jLiparill2006l ). The 
background induced by the atmospheric neutrinos is also shown in the figure, which is from a 
102 deg^ sky region with two 3° x3° areas around the two TeV sources, TeV J2021+4130 and 
MGRO_J2019+37, exclu ded (lAbdo et al.ll2007al ). The atmospheric neutrino flux is adopted 
from iHonda et al.l (120071 ) , with the zenith angle cos -0.6 -0.7 B 



It can be seen from Fig. [3] that if the cut-off energy of protons is high enough, the 
neutrino induced muons are marginally observable over the atmospheric background for 
energy > 10 TeV. For a 15-year observation of IceCube, we estimate the significance ~ 
3a for an energy threshold 10 Te V. Compared with the neutrino emission from MGRO 
J2019-I-37 (IBeacom fc Kistlerl 120071 ). our result of diffuse neutrinos is ~ 4 times larger. This 
is reasonable as the total diffuse TeV p hoton emission is about 4 times larger than that from 
MGRO J2019+37 (lAbdo et al.ll2007al ). However, since the diffuse atmospheric background 
is abo ut 30 times higher than the case of MGRO J201 9-I-37, the detectability becomes 
worse (IBeacom &: Kistlerl 120071 ) . iKistler fc BeacomI (120061 ) got the similar conclusion about 
the diffuse neutrinos. If the cut-off energy of protons is as low as ~ 1 PeV, the signal is 
dominated by the atmospheric background and is almost invisible. 

It should be pointed out that the neutrinos from the GCRs cosmic rays interacting with 
the ISM are negligible and not shown here. From Fig{T]we can see that the 7r° component 
of 7-ray flux from GCR background is about one order of magnitude lower than the Milagro 



•^The zenith angle of the Cygnus region relative to the detector at the South pole is about 130°, i.e., the 
neutrinos are up-going. Since there is a mirror symmetry between the up-going and down-going atmospheric 
neutrinos, it is equivalent to the atmospheric neutrino at direction with zenith angle ~ 50°. 



- 8 - 



data. Therefore the neutrino emission from this component is also very small. It can also 
be noted that the diffuse 7-ray emission is actually not uniformly distributed in the Cygnus 
region, but concentrates in some areas. Therefore the atmospheric background will be lower 
if proper observational region is concentrated, and then the detectability of signals will be 
improved. 




Fig. 3. — The cumulative muon events induced by neutrinos from the Cygnus region and 
atmospheric background on IceCube for one-year observation. 



4.2. The leptonic model 



A high energy electron population similar with the hadron ic one is also able to produce 
the TeV 7-rays through IC scattering (lAharonian et al.l l2006l ) . Similar with the hadronic 
model, an electron source population with differential energy spectrum d0e/d-Ee E~"^ exp{—Ee/E^) 
is introduced to account for the TeV 7-ray emission. The high energy spectral index of the 
electron population is adopted as ag = 2.6, which is harder than the Galactic background 
electron spectrum ~ 3.3. The cut-off energy is adopted as ~ 100 — 1000 TeV. Note that 
due to the fast energy losses of electrons in interstellar magnetic fi eld, the maxim um energy 
of electrons generally cannot be higher than several hundred TeV (iGaisserl Il990l ) . 



However, such an electron population will contribute too many low energy 7-rays 
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through bremsstrahlung radiation and will exceed the EGRET data. Therefore a broken 
power-law of spectral index is introduced to match the data. W e adopt Oip = .l-5 for energy 



lower than 4 GeV, which is similar to that done in GALPROP (IStrong et al.ll2004a( ) 



For the high energy slope, it is not arbitrary due to observations and theoretical ar- 
guments. It is known that the fast energy loss through synchrotron and IC will soften the 
power index of electron spectrum by about 1. Therefore if the source spectrum is ~ 2 as 
expected from 1st order Fermi acceleration, the propagated one is ~ 3. However, we argue 
that the reacceleration of electrons by the stellar wind shock of massive stars in the Cygnus 
region can maintain a harder spectrum — 2.6 as we adopted. 

The target ISRF field is adopt ed as the corn bination of three components: CMB, infrared 
from dust and optical from stars (jPorterll2005l ). The predicted 7-ray spectrum is shown in 
FiglH It shows that this model can also give good description to data. The total energy of 
this electron source is estimated as W{> 1 GeV) ~ 10 



49/ 



1 kpc ' 



erg. 




Fig. 4. — The energy spectrum of diffuse 7-rays in the Cygnus region as measured by 
EGRET and Milagro and predicted by the GALPROP model combined with additional 
electron source and DM annihilation components. See text for details. 



It is expected that such a high energy electron population will generate low energy 
synchrotron radiation in the interstellar magnetic field. We show the synchrotron radiation 
in Figl5]for the magnetic field -B = 1, 3 and 10 fiG respectively. The synchrotron spectrum 
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is ~ + l)/2 = 1.8, which is significantly different from the one by GCR background 
electrons, as shown in FigjSl 

Compared with the GCR background contribution, the additional electrons will con- 
tribute high energy synchrotron radiation, and might be detected by X-ray det ectors. We find 
in th e soft X-ray background (SXRB) from ROSAT All-Sky Survey (RASS, ISnowden et al. 
19971 ) that there is no distinct excess in the Cygnus region compared with nearby sky regions. 
The result from RASS in the Cygnus region is shown in this figure as an upper limit of the 
synchrotron radiation. It is shown that even for the magnetic field 10 /iC the predicted 
synchrotron radiation in the soft X-ray band is consistent with the observational SXRB. 
However, the peak of the synchrotron radiation is in the hard X-ray band (~ hundred keV) 
if the electron source has an energy cut ~ 1 Pe V. It may be detectable by the future satel- 
lite Hard X-ray Modulation Telescope (HXMT; Iwu et al.l[2002l : lli et al.lEoOfil ). HXMT is 
a space satellite with high sensitivity and a field of view 5.7° x 5.7°, which is devoted to 
performing a hard X-ray all-sky imaging survey. The wide field of view and high sensitivity 
make it possible to detect the diffuse X-ray emission from the Cygnus region. The sensitivity 
curve of HXMT for the fix-direction observation of diffuse emission is shown in FigjSl If the 
astrophysical background in hard X-ray band is not too high, the diffuse emission in the 
Cygnus region is detectable. If = 100 TeV, the "TeV excess" of Milagro data can also be 
reproduced, but the synchrotron radiation is almost invisible due to the high SXRB. 



5. Conclusions and Discussions 

The TeV 7-ray emission of the Cygnus region observed by Milagro shows significant 
excess compared with the conventional CR model. In this work we introduce a high energy 
proton or electron source to explain this high energy 7-ray emission. It is shown that both 
the hadronic and leptonic models can give good explanation to the current data. The total 
energy needed for the source population is consistent with the typical energy output of a 
supernova. 

The associated neutrino emission for the hadronic model and synchrotron radiation for 
the leptonic model are discussed as possible probes to discriminate these two scenarios. For 
the hadronic model, we estimate a 3a significance for 15-year exposure of IceCube in case 
the protons cutoff energy reaches ~ 10 PeV. The non-null result of the neutrino detection 
will support the hadronic model. For the leptonic scenario, the synchrotron radiation by the 
electrons is more luminous at the X-ray wavelength. We show that if the energy cut-off of 
the electron source is as high as ~ 1 PeV, the synchrotron radiation will peak in the hard 
X-ray bands and it might be observable on the satellite HXMT. 
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Fig. 5. — The synchrotron radiation from the electron source population for magnetic 
fields B = 1, 3 and 10/xG (from bottom to top) respectively. The GCR background syn- 
chrotron radiation is calculated using GALPROP. The observed X-ray limit is from RASS 
(ISnowden et al.l 119971 ). Also shown is the sensitivity curve for fix-direction observation of 
HXMT. 



In the calculation, the ISM density, the ISRF intensity and the magnetic field are 
adopted as typical Galactic values. Since there are a large number of massive stars and 
molecular clouds in the Cygnus region, the ISRF and ISM density may be higher than the 
Galactic average values. We should point out that this does not change the basic conclusion 
in this work. This is because a higher ISM density or ISRF intensity can be effectively 
compensated by a lower CR fiux. For the hadronic model, the neutrino fiux will keep the 
same since neutrino fiux is proportional to the 7-ray fiux. For the leptonic model a higher 
ISRF intensity will correspond to a lower electron luminosity. However, considering the 
magnetic field in the Cygnus region may be also stronger than the Galactic value we may 
expect similar synchrotron radiation as given in the work. 

Finally, since the present data are lack in the energy interval from ~ lOGeV to ~ 
ITeV due to the transition of detection technology from space to ground there are large 
uncertainties for model construction. The space telescope F ermi will extend the detection 



energy up to ~ SOOGeV (IRitz fc GLAST Mission Teamll2004j ). Together with some ground- 



based experiments , such as ARGO, which may lower the energy threshold down to ~ lOOGeV 
(IBacci et al.lll999l ). we will have better understanding in the emission mechanism of this 
region in the future. 



We thank Lu Ye, Lu Fang-Jun, Qu Jin-Lu and Zhang Shuang-Nan for helpful discus- 
sions. The work is supported by the NSF of China under the grant Nos. 10575111, 10773011 
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